Generalized Phase Field Model for Computer Simulation of Grain Growth in 

Anisotropic Systems 



0^ 
o 

Q 

(N 



C/2 



o 



> 
o 

On 
ON 



I 

O 

o 



A. Kazaryan^, Y. Wang^, S.A. Dregia^ and Bruce R. Patton^ 

^Department of Physics, ^Department of Materials Science and Engineering, The Ohio State University, Columbus, OH 43210 

(February 1, 2008) 

We study the dynamics and morphology of grain growth with anisotropic energy and mobility 
of grain boundaries using a generalized phase field model. In contrast to previous studies, both 
inclination and misorientation of the boundaries are considered. The model is first validated against 
exact analytical solutions for the classical problem of an island grain embedded in an infinite matrix. 
It is found that grain boundary energy anisotropy has a much stronger effect on grain shape than 
that of mobility anisotropy. In a polycrystalline system with mobility anisotropy, we find that the 
system evolves in a non-self-similar manner and grain shape anisotropy develops. However, the 
average area of the grains grows linearly with time, as in an isotropic system. 

PACS numbers; 68.55.-a, Sl.lSAa, 05.45.-a 



The effect of anisotropy in energy and mobility of grain 
boundaries on the kinetics of grain growth and mor- 
phological evolution is a relatively unexplored problem. 
Studies have shown that energy and mobility of a grain 
boundary depend on the misorientation between the two 
crystals and the inclination of the grain boundary In 
addition, phenomena as segregation of impurities l2] or 
presence of a liquid phase ^ at the grain boundaries may 
also result in anisotropy of both energy and mobility. 

Although most computer simulations of grain growth 
have been performed for isotropic cases , anisotropy 
in grain boundary properties has been introduced in 
a number of simulations, mainly by the Monte Carlo 
method p|-p^. For example in the study of texture de- 
velopment during grain gr owth grains are divided 
into two types and the contacts between them form three 
kinds of grain boundaries of either small or large misori- 
entations. To take into account the full range of grain 
orientations, more general approaches have been pro- 
posed Ip^-p^. However, all these models consider either 
misorientation or inclination dependence of grain bound- 
ary properties. A simple dislocation model of the grain 
boundary shows that both energy and mobility of the 
boundary could depend strongly on both misorientation 
and inclination fl^ . Furthermore, in the models that 
deal with inclination dependence only a few incli- 
nations are considered and the use of only first neighbors 
for the calculation of boundary inclination results in an 
intrinsic energy and mobility anisotropy associated with 
the discrete lattice used in the simulations . 

In this paper, to study the kinetics and morphology 
of grain growth in anisotropic systems, we extend the 
phase field approach to take into account both inclina- 
tion and misorientation dependence of grain boundary 
energy and mobility. The phase field method has been 
successfully applied for computer simulation of isotropic 
grain growt h ||7| ,p[|, phase transformations and so- 
lidification ]19||. In this model, the polycrystalline mi- 



crostructure is described by a set of non-conserved order 
parameter fields (?7i, 7^2, ... , r]p), each representing grains 
of a given crystallographic orientation. Microstructural 
evolution of a polycrystalline system is characterized by 
the spatio-temporal evolution of the order parameters, 
via the Ginzburg-Landau-type kinetic equations: 
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dt 
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where L is the kinetic coefficient characterizing the grain 
boundary mobility, and F is the free energy functional of 
the form: 
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where fijii) is the local free energy density and k is the 
gradient coefficient, which together determine the width 
and energy of the grain boundary regions. 

Several attempts have been made to introduce 
anisotropic boundary properties into the phase field for- 
mulation. For example, in the phase field models of 
antiphase domain growth and solidification, anisotropy 
has been introduced by using multiple order parameters 
based on the underlying crystal symmetry po| , pT| . How- 
ever, application of these models to grain growth is very 
difficult due to the fact that the description of the grain 
boundary is much more complicated, and a quantitative 
description of a general grain boundary is still lacking. 
Recently, Kobayashi et al. ||2^ and independently Luck 
psf suggested a different way to extend the phase field 
model to include anisotropy in grain boundary proper- 
ties. In their approach, anisotropy is incorporated by 
introducing an additional variable that describes spatial 
orientation of the grains. So far, the model has been 
applied to only quasi-one-dimensional systems with mis- 
orientation dependence of grain boundary energy. 
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In contrast, a simple phenomenological approach has 
been successfully used to describe the surface energy and 
mobility anisotropy of crystal surfaces in phase field mod- 
eling of crystal growth, where the gradient coefficient k 
and the kinetic coefficient L have been formulated as 
functions of crystal surface orientations jl^j2^. In the 
current paper a similar approach is used to describe grain 
boundary energy and mobility anisotropy. For the sake 
of simplicity, we consider tilt grain boundaries between 
two-dimensional crystals with square lattices, which (ig- 
noring rigid-body translations) can be described by two 
parameters: relative orientation 9 of the grains (misori- 
entation) and spatial orientation of the boundary with 
respect to the reference coordinate system (inclination). 
For crystals with four-fold symmetry and based on sim- 
ple dislocation models |]l6| , the energy of a low-angle tilt 
boundary {9 < 20") can be approximated as: 

E{9, 4,) = Eoi\cos{(j))\ + \sini4)\)9{l - ln{9/9^)) (3) 

where = corresponds to the symmetric tilt boundary 
(see Fig.^, 9m is the misorientation at which energy is 
maximum and Eq is a constant. The grain boundary 
energy can be related to the gradient coefficient by k (x 
E'^{9, (j)) However, as noted by McFadden et al, the 
gradient coefficient must be a differentiable function with 
respect to inclination for the Ginzburg-Landau equations 
to be properly defined p^ . Thus in our simulations the 
grain boundary energy has been taken in the form: 

E{9, 4,) = Eo{l - fecos(40))0(l - ln{9/9m)) (4) 

which maintains four-fold symmetry in 0, where 6e is 
a phenomenological parameter serving as a measure of 
the degree of anisotropy. The inclination dependence of 
the energy is similar to the one used in the phase field 
model of solidification ||l9|. By choosing 6e — 0.24 the 
same anisotropy ratio as in Eq. (^) is obtained, where 
the anisotropy ratio r is defined as the ratio of the largest 
to the smallest grain boundary energy at a fixed misori- 
entation ^o- f = maxE((l),9Q)/minE{(j),9Q). 

In the phase field approach to grain growth, grain 
boundary mobility is characterized by the kinetic coef- 
ficient L. In contrast to the grain boundary energy, val- 
ues of the grain boundary mobility are harder to esti- 
mate. Qualitatively, in pure materials the higher the 
defect concentration at a grain boundary the higher its 
energy as well as its mobility. To investigate qualitatively 
the effect of mobility anisotropy on the behavior of grain 
growth, we assume mobility to have the same dependence 
on misorientation and inclination as the grain boundary 
energy (Eq.(^), with a different phenomenological pa- 
rameter Sl. By varying Se and Sl we can control the 
anisotropy ratio and investigate the interplay between 
energy and mobility anisotropy. 

The free energy density in Eq. (^ has been taken in the 
form 0,|: 



i—l 2—1 j>i 

In our simulations grain boundary misorientation was 
considered to be in the range of < 6* < 0, where 
Q < 9m Then, misorientation 9ij between grain i and 
grain j = I . . .p, where p is the total number of or- 
der parameters used in the simulation) is calculated as 
\i — j\ ■ Q/{p — 1). The following values for the phe- 
nomenological parameters in Eq.(|^) have been used in 
the simulations: Oi = 1.0, 02 = 1.0, 03 = 2.0. Eq.([T]) was 
discrctizcd using the second order Euler technique on a 
unit square lattice with dt = 0.1. 

Below we present several applications of our model to 
systems containing an island grain in an infinite matrix 
as well as to a polycrystalline aggregate. 

Island grain with energy anisotropy. To validate the 
model against exact analytical results, we first examine 
shrinkage of an island grain embedded in an infinite ma- 
trix. If the grain shrinks in a "self-similar" manner, i.e. 
maintains its shape while shrinking, Taylor and Cahn ]2^ ] 
have shown that the Ginzburg-Landau equations could 
be rewritten in the same form as in the isotropic case. As 
a result, shrinkage kinetics of a grain with arbitrary but 
self-similar shape should be the same as in the isotropic 
case. Moreover, they have argued that depending on 
the form of the gradient coefficient certain complications 
could occur in the solution of the Ginzburg-Landau equa- 
tions due to non-convexity in the Wulff plots. Since 
we are mainly interested in the qualitative behavior of 
the system and to avoid those complications, we have 
chosen a grain boundary energy in the form (^ with 
Se = 0.05 and Eq = 1/(1 + Se), which gives anisotropy 
ratio rE — 1.1. The result of our simulations is presented 
in Fig.^. As expected, an initially circular grain trans- 
forms to its Wulff shape and then shrinks in a self-similar 
manner. The extracted kinetics indeed shows linear de- 
cline of area with time, in agreement with the theoretical 
predictions of Taylor and Cahn |2^] . 

Island grain with mobility anisotropy. Following the 
theoretical analysis of Allen and Cahn ||2^ for the 
isotropic case, we obtain a similar relationship for the 
shrinkage rate of a single island grain with anisotropic 
boundary mobility: 

A{t) = A{0) - 2TTk{L)^t (6) 

where A is the area of the grain, and (. . .)^ represents 
an average over all possible inclinations. This expres- 
sion is valid for arbitrary particle shapes and the sys- 
tem does not need to be in a self-similar regime. In 
the simulations mobility has been taken in the form (Q) 
with Lq = + Sl). In contrast to the case of en- 
ergy anisotropy, a small mobility anisotropy ((5l=0.05, 
which gives an anisotropy ratio Lmax/Lmin = 1-1) has 
little effect on the shapes of the grain. Indeed, to reach 
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a similar degree of anisotropy in shape(Fig.^), we have 
to introduce a significantly higher mobility anisotropy 
with 6l—0-7, which produces a relative mobility differ- 
ence {{Ljnax - Lmin) / Lmin) ~ 500% (while Only 10% dif- 
ference in energy is required) . Experimental observations 
have shown that the mobility difference could be orders of 
magnitude larger than the energy diffrence for the grain 
boundaries [p]j9|JlO[| . Finally, excellent agreement between 
the simulation results and the exact analytical solution 
(|) for the grain shrinkage rate has been observed (FigJ^), 
which indicates that discretization of the equations does 
not influence the kinetics. 

Island grain with both energy and mobility anisotropy. 
Next, we discuss the effect of the interplay between the 
energy and mobility anisotropy of the grain boundary on 
the shape of a single shrinking island grain. It should be 
noted that since the functions E{(j)) and L{(p) have been 
taken to have the same (j) dependence, for a shrinking 
island grain energy and mobility anisotropy have oppo- 
site effects on the shape. Boundaries with normals in 
■|n directions have the lowest energies as well as lowest 
mobilities. Therefore, energy minimization would prefer 
a shape bounded by the low energy boundaries (Fig|2^) , 
while grain shrinkage with anisotropic boundary mobility 
would produce a shape bounded by the fastest moving 
boundaries (Fig.^lb). As we mentioned earlier, a small 
anisotropy in the mobility coefficient {r^ = 1-1) does not 
have much effect on the shape of the grain; thus the grain 
shape is dominated by energy anisotropy, giving a shape 
similar to Figj^a. On the other extreme a very large mo- 
bility anisotropy (r^ = 40) results in diamond-like shapes 
(Fig.^). In the intermediate range (r£;=l.l, ri=5.7), a 
grain shape with quasi eight-fold symmetry has been ob- 
served (Fig.||:), although both energy and mobility as a 
function of inclination have four-fold symmetry. 

Poly crystalline aggregate with mobility anisotopy. For 
a polycrystalline aggregate, MuUins has shown that the 
average grain area grows linearly with time if the system 
evolves in a statistically self-similar (scaling) manner (i.e. 
all configurations have identical statistics when trans- 
formed to the same linear scale by uniform magnification) 
p7| . However, self-similarity may not hold when bound- 
ary properties (energy and/or mobility) are anisotropic 
p8| . Indeed, our simulations performed on a polycrys- 
talline system with grain boundary mobility anisotropy 
have shown that grain shapes evolve in a non-self-similar 
manner, e.g., shape anisotropy develops (Fig.^). The 
simulations were performed on 512 x 512 square lat- 
tice with 36 order parameters. The simulations were 
started from an isotropic polycrystalline microstructure 
consisting of ~ 1000 grains, obtained from nuclcation and 
growth of crystals from a liquid phase in an isotropic sys- 
tem. Grain boundary mobility anisotropy is introduced 
according to L = 9{1 — ln{9/dm)){l — dcos{2(j))), with 
0^ = 10°, 9 = 5° and Sl = 0.9, which has two-fold 
symmetry with respect to </> and (j) = n/2 is the fastest 



growth direction. For comparison, the result obtained 
from the same initial microstructure in an isotropic sys- 
tem is shown in Fig.|^. Comparing these results it is clear 
that in the anisotropic case grain shape anisotropy devel- 
ops. Quantitative analysis of Fig.^ has shown that the 
inclination distribution of grain boundaries (i.e., bound- 
ary length as a function of inclination) is no longer uni- 
form, which suggests that the microstructure is no longer 
self-similar. However, the average grain area still grows 
linearly with time as shown in Fig.|| It can be shown an- 
alytically that even though self-similarity is not sat- 
isfied when grain boundary mobility is anisotropic, the 
average grain area may still grow linearly with time. 

In conclusion, we have studied the dynamics and mor- 
phology of grain growth with both energy and mobility 
anisotropy of grain boundaries using a generalized phase 
field model that incorporates both inclination and misori- 
entation dependence of grain boundary properties. For 
a single island grain embedded in an infinite matrix we 
found that a small anisotropy in grain boundary energy 
can produce strongly anisotropic grain shapes. In con- 
trast, the mobility anisotropy needs to be significantly 
stronger to produce similar morphologies. For a poly- 
crystalline aggregate with mobility anisotropy, we found 
that the average area of the grains grows linearly with 
time, as in the isotropic case, even though grain shape 
anisotropy develops, which indicates that the microstruc- 
tural evolution is no longer self-similar. 
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FIG. 1. Schematic drawing of a tilt grain boundary formed 
by two crystals with misorientation Q and inclination (\), mea- 
sured from the symmetric tilt boundary 



(c) 

FIG. 2. Shrinkage of a single island grain with (a) energy 
anisotropy; (b) mobility anisotropy; (c) both energy and mo- 
bility anisotropy 
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FIG. 3. Shrinkage kinetics of an island grain embedded in 
an infinite matrix; comparison between analytical solution 
Eq.(^ (solid curve) and simulation (solid circles) in the case 
of anisotropic grain boundary mobility. Ao is the initial area 
of the grain. 




(a) (h) 
FIG. 4. Typical microstructures obtained from computer 
simulations with isotropic grain boundary energy and (a) 
isotropic and (b) anisotropic boundary mobility with two-fold 
symmetry with respect to 
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FIG. 5. Time dependence of the average grain area 
in a polycrystalline system with grain boundary mobility 
anisotropy. 



4 



